Research in Astron. Astrophys. 20xx Vol. 9 No. XX, 000-000 

http://www. ma-journal, org http://www. iop. org/journals/raa 



Research in 
Astronomy and 
Astrophysics 



Detailed Analysis of Fan-Shaped jets in Three Dimensional 
Numerical Simulation 

R. L. K. Shibata^, H. Isobe^ and C. Fang^'^ 

^ Department of Astronomy, Nanjing University, Nanjing 210093, China; rljiang@nju.edu.cn 
^ Kwasan and Hida Observatories, Kyoto University, Yamashina, Kyoto 607-8471, Japan 
^ Unit of synergetic Studies for Space, Kyoto University, Yamashina, Kyoto 607-8471, Japan 
Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of 
Education, China Received [year] [month] [day]; accepted [year] [month] [day] 



Abstract We performed three dimensional resistive magnetohydrodynamic simulations 
to study the magnetic reconnection using an initially shearing magnetic field configura- 
tion (force free field with a current sheet in the middle of the computational box). It is 
shown that there are two types of reconnection jets: the ordinary reconnection jets and 
fan-shaped jets, which are formed along the guide magnetic field. The fan-shaped jets are 
much different from the ordinary reconnection jets which are ejected by magnetic tension 
force. There are two driving forces for accelerating the fan-shaped jets. The one is the 
Lorentz force which dominates the motion of fluid elements at first and then the gas pres- 
sure gradient force accelerates the fluid elements in the later stage. The dependence on 
magnetic reconnection angle and resistivity value has also been studied. The formation 
and evolution of these jets provide a new understanding of dynamic magnetohydrody- 
namic jets. 

Key words: magnetohydrodynamics: MHD — methods; numerical 
1 INTRODUCTION 

Magnetic reconnection plays a very important role in solar flares, coronal mass ejections, and other solar 
activities. The problem of flare energy release has been puzzling people for many years before the devel- 
op ment of magnetic r econnection m echanism. The first pioneers in this field of magnetic reconnection 
are lSwee3 (Il958b and lParkeJ (Il957l) . who developed a theory well known as Sweet-Parker mechanism. 
However, it can not explain the solar flares. The biggest problem in their theory is that the time scale 
of energy release is much longer than the realistic one of flares. Several years later, [Petschek (1964j) 
improved this work by considering a pair of slow mode shocks outside a very small diffusion region, 
which greatly increases the reconnection rate. The fast reconnection model is successful, though some 
people do not believe that the large energy release p rocess has such a small diffusion reg ion. The famous 
CSHKP model dCarmichaelll 19641; ISturrock|[T966t lHiravama|[l974l; iKopp & Pneuma nl ll976h based on 
magnetic reconnection for flares has been developed since 1960s. 

Besides observatio ns, magnetohydrodynamic (MHD) simulations play a key role in studying the 
magnetic reconnectio n (Ugai & Tsuda 1977; Forbes & Priest 1983; Forbes 1990; iForbes & Malherbel 
ll99lHYokovama & Shibata 1994; Ma gara et al. 1996; Ugai 1996; Chen et al. ,1999D. With flie develop- 
ment of computer and observation technology, the CSHKP flare model has been extended. More and 
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more simulations show that the plasmoid ejection may greatly increase the mag netic reconnection rate 
dShibata et al.ll 19951: IShibatall 19961 1 19971 iMagara et al.lll997l; iNishida et a[mO% . The plasmoids can 
form in an anti-parallel magnetic field. Before ejection magnetic energy is restored in the diffusion re- 
gion, and simultaneously the plasmoids merge with each other Finally, the restored magnetic energy is 
released in a very short time scale after the ejection of the big plasmoid. 

Recently, the three dimensional (3D) simulation is one of the hottest points in astrophysics. In this 
case, the reconnecting process is much more complicated than th at in two dimensional (2D) case. The 3D 
simulations can provide us very complex and realistic structures ([Yokovama & Shibatall995lllsobe et al.l 
l2005l lUgai & Zhen"gll2005t llsobe et al.ll2006 l: 'Shimizu et al."2009V It is also a good tool to explain the 
mech anism of many small scale solar activities (Cirtain et al. 2007; Shibata et al. 2007; Katsukawa et af] 
I2OO7I) which have been observed by Hinode dKosugi etal.r2007h . Some of these observed jetlike features 
are believed to be a result of magnetic reconnection in the so-called interlocking-comb like magnetic 
configuration in sunsport penumbra (i.e., penumbral jets). The motivation of this paper is to study the 
effect of the guide field using 3D numerical simulations. In our new simulation, we found that some 
parts of the reconnection jets which are ejected from the diffusion region can move along the magnetic 
guide field lines in the interlocking-comb like magnetic configuration. Moreover, the reconnection jets 
which move along the guide field lines are almost perpendicular to the ambient magnetic field. In order 
to distinguish the jets moving in different directions, we call the jets moving along the ambient magnetic 
field as ordinary reconnection jets and the jets moving along the guide field as fan-shaped jets (we will 
see the shape of the jets is similar to a fan in the next section). Besides our simulations, two very recent 
paper by Pontin et al. (2005) and Ugai (2010) also discussed the effect of guide field. Their simulations 
showed that the guide field can distort the propagation of 3D plasmoids and the plasma is almost ejected 
along the guide field. It seems that the guide fiel d plays a very im portant role for forming these jets in 
the 3D magnetic reconnection. In our first paper (jjiang et al.ll201 ll) we briefly reported the main results 
of our 3D simulation, while in this paper the detailed analysis and discussion of fan-shaped jets are 
presented. 

This paper is organized as follows. A detailed description of the magnetohydrodynamics (MHD) 
equations and the initial condition are given in the next sections. In Section[3]we show the 3D, 2D and 
ID structures of our typical simulation results. Moreover, the analysis of the Lagrangian fluid elements, 
reconnection rate and slow mode shock are also given. Section|4]shows dependence on different recon- 
nection angles (defined in the initial conditon) and different values of resistivity in the diffusion region. 
Finally, discussion and summary are given in the last Section. 

2 NUMERICAL METHOD 

We solve the three-dimensional, nonlinear, resistive, compressible MHD equations numerically. The role 
of gravity, thermal conduction, radiation or some other effects are not so important in our simulations, 
because what we are really interested in is the dynamics of the magnetic reconnection. The additional 
consideration may affect the final results but not essential. The simplified MHD equations are given in 
cgs units as follows: 
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where the eight independent variables are the velocity {vx,Vy,Vz), density (p), magnetic field 
{Bx, By, Bz), and pressure {p). rj is the resistivity, j the current density, T the temperature, and E the 
electric field. The specific heat ratio of monoatomic perfect gas is 7 = 5/3. The normalized units of the 
variables are listed in Table [T] Hereafter, all the values used in our analysis and discussion are in non 
dimensional form for simplicity. 



Table 1 Normalization Units. 



Variable 



Quantity 



Unit 



x,y,z 


Length 


u 


T 


Temperature 


To 


P 


Density 


po 


P 


Pressure 


Po = poTaKB/m 


t 


Time 


to = Lo/vo 
770 = {l/c^)LoVo 


V 


Resistivity 


V 


Velocity 


vo = (polpoY''^ 


B 


Magnetic field 


D 1/2 
Bo = Po 


E 


Electric field 


Eo = {l/c)voBo 


j 


Current density 


jo = {c)Bo/Lo 


F 


Force per unit mass 


Fa = po/Lo 


VA 


Alfven velocity 


VAO = Bo/{poY^'' 


TA 


Alfven time scale 


TAO = Lo/VA 


Ea 


Alfven electric field 


Eao = [1/c)vaBo 



Notes: In this table, m, kb and c mean molecular mass, Boltzmann constant and light speed, respectively. 
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Reconnection angle 



Fig. 1 Initial condition, and the box size is —4 < a; < 4, —16 < y < 16, — 16 < z < 16 in 
non dimensional unit. 



The computational box is shown in Figure[T] The domain of this box is —4 < .t < 4, — 16 < ?/ < 16 
and — 16 < z < 16 in all of our cases (Note that the images of simulation results only show the central 
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part of the computational box and the scale for x, y and z coordinates is different because of the thin 
current sheet in x direction). The most important condition is the configuration of the force free magnetic 
field which is shown in the right panel of Figure [T| The mathematical form of this field is 

= 0, (8) 



-Bmt fox X < -Ahx 

By = { B„„ sin {[0a; - (tt - 0) A/i,]/2A/i,} for \x\ < Ah, , (9) 
Bini sin [{20 - tt) /2] for x > Ah^ 

for X < —Ahx 

B, = { B„„ cos {[0a; - - e)AK]/2Ah,} for |x| < Ah, , (10) 
Bim COS [{29 - 7r)/2] for x > Ah, 

where Bini is a constant, which means the magnetic field has an initial uniform strength. The direction 
of this field is shown in Figure[T] The Ah, is the half width of the current sheet. 9 is the reconnection 
angle shown in Figure[T] and it has been taken to be one of the four values in our simulation cases, i.e., 
TT, 37r/4, 7r/2 and 7r/3, as listed in Table |2l 



Table 2 Computational Cases. 


Case 


/3 




Reynolds number (vaLx /rjini) 


Reconnection angle (9) 


Total grid points 


1 


1.0 


0.025 


226 


TT 


300 X 240 X 240 


2 


0.8 


0.050 


126 


TT 


240 X 480 X 480 


3 


0.8 


0.050 


126 


37r/4 


240 X 480 X 480 


4 


0.8 


0.050 


126 


tt/2 


240 X 480 X 480 


5 


0.8 


0.050 


126 


n/3 


240 X 480 X 480 


6 


1.0 


0.025 


226 


n 


240 X 480 X 480 


7 


1.0 


0.050 


113 


TT 


240 X 480 X 480 


8 


1.0 


0.075 


75 


TT 


240 X 480 X 480 



Notes: In this table, Lx is the half length of the computational box size in x direction. 



In order to get fast reconnection process, we put a localized diffusion region at the center of the 
current sheet ( fUg ai 1986; ,Yokoyama & Shibata,1994) . The adopted functional form for the anomalous 
resistivity is as follows: 



XTT XIT 

T] = '^^^ cos[^^^] cos[^— -] cos[^^^] for |a;| < Ah,, [y[ < Ahy, [z[ < Ah^, (11) 

X y z 

where Ah, = 0.4 (it is the same value as the width of the current sheet), Ahy = 1.6, Ah^ = 1.6 is 
the half width of the diffusion region in x, y and z direction, respectively. 77,„, is the maximum of the 
resistivity in the diffusion region and is taken as a constant. 

In the numerical procedures, we solve the MH D equations using CIP- MOCCT scheme 
dYabe & Aokilll99ll:lYabe et alJll99ll:lKudoh et al.lll999h with an artificial viscosity dStone & NormanI 



119921) and libraries of CANS (Coordinated Astronomical Numerical Softwares). The calculation do- 
main is discretized into 240 x 480 x 480 nonuniform grids in most of cases. The minimum grid size 
is Aa; = 0.00667, Ay = 0.0667 and Az = 0.0667. Since there is no gravity in our simulations, the 
free boundary condition is simply applied in all of six boundaries (namely, an equivalent extrapolation 
was applied for all quantities). The non- divergence condit i on of the magnetic field is guaranteed by 
the constraint transport (CT) algorithm of lEvans & HawlevI (Il988b in CIP-MOCCT scheme. For initial 
condition, the density and pressure are uniform (p = 1, p = 1) in the computational box and the other 
parameters are shown in Table |2] 
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3 TYPICAL CASE 

Our typical case is the one in which the reconnection angle is tt with a plasma beta /3 = 0.8 and 
Vini — 0.050 (case 2 in Table |2]i. Figure |2] shows the 3D visualization of gas pressure distribution at 
different times. The bright colors and high opacity stand for the high values. The volume rendering 
in the 3D visualization does not show the low gas pressure volumes which exist everywhere in the 
computational box. Magnetic field lines are drawn from Lagrangian fluid elements located outside of 
the diffusion region, and the same field lines are shown in the upper panels of Figure |2] We can see a 
fan-shaped structure in the current sheet. For the lower panels, no magnetic field lines are plotted but 
we add a cross section near z ~ and a velocity field plane at a; = 0. The velocity arrows show the 
fan-shaped structure almost along the guide field lines (there is no legend for the velocity field, the scale 
for the velocity is useless in these 3D images because of the projection effect). The cross section located 
near z = 0.0 shows the ordinary magnetic reconnection jets in t his plane which is very similar to the 2D 
or 2.5D (two dimensio n with three components) simulations jChen et al. 1 1200 It lYokovama & Shibatal 
I2OOII: IJiang e"taDl2010h . We will discuss this in detail in the next paragraph. Figure [3] shows the 3D 
results at different view points at the time = 11.5. From this figure we can see the fan-shaped structure 
jets more clearly. 

The 2D cross sections at different positions and different times are shown in Figures |4]|6l Figure |4] 
depicts the gas pressure, magnetic pressure and total pressure distributions inx — y plane at z = 0.0 (we 
can also briefly understand the position by the small bos above the color bar in this figure). The guide 
field is perpendicular to this plane, which means what we present in this plane is similar to simulations 
in 2D or 2.5D case as we mentioned above. As reconnection goes on, the central diffusion region is 
heated by Joule dissipation, and a pair of ordinary reconnection jets are ejected by magnetic tension 
force. The central low magnetic pressure region is due to the magnetic reconnection process, which can 
be compared to the high gas pressure region. For a steady state, the two pressures are almost balanced 
with each other in the x direction of the inflow region shown by the total pressure in the lower panels 
of Figure |4l The X-shaped structures (we can find them in the panels showing the gas and magnetic 
pressures) almost disappear. 

Figure|5]shows the gas pressure, magnetic pressure and total pressure distributions in y — z plane at 
a; = 0.0. In this plane, the ordinary reconnection jets are almost perpendicular to the magnetic field lines 
(in other words, moving along the y and —y directions) and the jets almost moving along the magnetic 
field lines (along the z and — z directions) are shown only in a cross-section of the fan-shaped jets. The 
velocity of the fan-shaped jets is smaller than that of the ordinary reconnection jets whose velocity is 
approximately the Alfven speed (about 1.58). High magnetic pressure regions have also been ejected as 
shown in the magnetic pressure distribution images in Figures |4] and |5] (high magnetic pressure parts, 
yellow color). The temperature and density distributions in this plane are shown in Figure |6l Similar to 
the gas and magnetic pressure distributions, the distributions of the density and temperature are roughly 
opposite with each other 

Figures |7] and [8] depict pressure distributions in y — z plane scanning at x = —0.1, x — 0.0 and 
X = 0.1 at the time = 8.5 and 11.5, respectively. Because the domains of these figures are in the current 
sheet, we can see that the magnetic field is shearing between different columns. The interesting thing is 
that the jets also have different directions for different x positions. In other words, the direction of the 
jets continuously changes following the direction of the magnetid field which is shearing in the current 
sheet for forming fan-shaped jets as shown in the 3D visualization in Figures |2] and |3] The cartoon in 
Figure |9] shows both ordinary reconnection jets and the fan-shaped jets. One can clearly understand the 
shape and the stucture of the jets by this cartoon. 

The interesting ID velocity distributions along x (y ~ 0.0, z = 0.0), y (x = 0.0, z = 0.0) and z 
{x = 0.0, y = 0.0) directions at time = 8.5 and 11.5 are shown in Figure[TO] Note that the velocities 
in X, y and z directions stand for inflow speed, ordinary reconnection outflow speed and fan-shaped 
jet speed, respectively. The inflow speed [x direction) becomes a little smaller at the later stage (time 
= 11.5), which indicates that the magnetic reconnection rate (as discussed later) reaches the maximum 
around the time = 8.5. At the latter stage the ordinary reconnection outflow speed is approximately the 
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Time=8.50 Time=11.5 




Fig. 2 3D visualization of the gas pressure distributions at different times (typical case). The 
bright colors and high opacity stand for the high gas pressure and the low value is transparent. 
Solid tubes are magnetic field lines in the upper panels and arrows are velocity in the lower 
panels. Since the current sheet is very thin, a different scale in x direction is used in these 
panels. 



ambient Alfven speed (1.58). Moreover, the maximum velocity of the fan-shaped jets along z direction 
is about 0.6. Thus, the driving forces of the fan-shaped jets and the ordinary reconnection outflow are 
different. 

In order to investigate the driving force of the fan-shaped jets, we examine the trajectories and the 
forces of Lagrangian fluid elements (test particl es). This is also a method to show how pla sma moves 
and what kind of force drives plasma to move dShibata & Uchidalll986l: [Kato et al.ll2002h . Figure [TT] 
shows the trajectory of a typical fluid element in the fan-shaped jet. This fluid element is located at the 
edge of the diffusion region at the initial time and the magnetic field is perpendicular to the trajectory 
(i.e., perpendicular to the velocity direction). After that, this element is ejected as a jet at the time = 8.5 
and eventually moves along the guide field line. The cosine value of the angle between velocity (v) and 
magnetic field (B) vectors is shown in the left panel of Figure[T2l At the initial time, the velocity of the 
particle is perpendicular to the magnetic field and almost parallel to that at the time = 14. The right 
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Time=11.5 Titne=11.5 




Fig. 3 3D visualization of the gas pressure distributions at different view points (typical case). 
The bright colors and high opacity stand for the high gas pressure and the low value is trans- 
parent. Solid tubes are magnetic field lines. 



panel of Figure [12] shows the forces in z direction on this element. It shows that there are two stages 
for acceleration. In the first stage (time = 8 — 10), the magnetic tension force drives this element since 
the gas pressure gradient force balances the magnetic pressure gradient force in the inflow region as 
we mentioned before. After that it is dominated by the gas pressure gradient in the later stage (time 
= 10 — 11.5). The trajectory in the x — z plane is shown in Figure [13] Note that the y direction is 
different at different times to follow the trajectory. At the time = 8.5, the element is driven as a inflow 
and at the time = 11.5 as a fan-shaped outflow. 

Because of the complexity in measuring the reconnection rate in 3D simulations, we give two kinds 
of measurement for the reconnection rate using the inflow speed and electric field as shown in Figure[T4l 
The inflow speed is the maximum one in x direction (at the positions y = and z = 0) and the electric 
field is taken at the point near the edge of the diffusion region {E^ and Ey are almost zero, so we can 
not see them in this figure). By comparing these two, one can have a overall view about how fast the 
magnetic field dissipates in our simulations. The upper panels in Figure [14] show that these two kinds 
of reconnection rate have a similar profile. Both of them reach the maximum almost at the same time 




Fig. 4 2D distributions (typical case) of the gas pressure, magnetic pressure and total pressure 
in a; — y plane at z = 0.0. The arrows in the upper panels indicate the magnetic field and others 
the velocity field. 
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Fig. 5 2D distributions (typical case) of gas pressure, magnetic pressure and total pressure in 
y — z plane at x = 0.0. The arrows in the upper panels indicate the magnetic field and others 
the velocity field. 
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Fig. 6 2D distributions (typical case) of density and temperature iny — z plane at a; = 0.0. 
The arrows in the upper panels indicate the magnetic field and others the velocity field. 



(around 8.5) as we analyzed in ID distributions and that is also the reason why we choose the time =8.5 
and a later time = 11.5 to show the 3D, 2D and ID distributions. The lower left panel in Figure fl4lis 
the mass and energy conservations with considering the flux coming in or out of the computational box. 
Even if CIP-MOCCT scheme is non-conservative for variables, our numerical error is still small enough 
to get a realistic result. Wt, Wk and Wm in the lower right panel are the total thermal energy, total 
kinetic energy and total magnetic energy, respectively. The magnetic energy release rate becomes faster 
and faster after the re connection rate re aching the maximum because of the formation and extension of 
the slow mode shockj Nitta et al.l200lh . Actually, most of magnetic energy is released by the slow mode 
shock CPetschek 1963^ 

The upper panels of Figure[T5]depict the gas pressure and current density distributions in z direction, 
which is only the central part of Figure |4] Other panels show one dimensional plot along the white line 
indicated in the upper left panel at the time = 14.0. We can see a pair of slow mode shocks in the 
current density distribution and the ID panels. The value jumps between two sides of the shock are 
shown by the shadow parts in the ID distributions. Maybe the shocks are not sharp enough, which is 
due to the lack of grid points (about 65 grid points along the white line). However, we can compare 
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Fig. 7 2D scanning figures (typical case) of the gas pressure, magnetic pressure and total 
pressure my — z plane at x = —0.1, 0.0, 0.1 and at time = 8.501. The arrows in the upper 
panels indicate the magnetic field and others the velocity field. 



this with the distribution of current density. The clear X-shaped structure is seen around the X-point. 
Moreover, the angle marked by a is related to the reconnection rate, bec ause we have t he formula 
tan(Q;/2) = Vin/vA according to the Petschek reconnection mechanism (|Petscheklll964l) . From the 
upper right panel, tan(a/2) is about 0.026, while the reconnection rate is approximate 0.032 given by 
the upper left panel of Figure[T4]at the time = 14.0. The reconnection rate is a bit larger than the value 
measured from the angle between two slow mode shocks. Considering the error of the measurement of 
the angle and the complexity in 3D MHD reconnection, we think that the two values are comparible. 
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Fig. 8 2D scanning figures (typical case) of the gas pressure, magnetic pressure and total 
pressure my — z plane at x = —0.1, 0.0, 0.1 and at time = 11.50. The arrows in the upper 
panels indicate the magnetic field and others the velocity field. 



4 PARAMETER DEPENDENCE 

In this section, we have two parameters to analyze. One is the reconnection angle (i.e., 6, shown in 
Figure [U, the other is the resistivity (i.e., ijini). Dependence on the reconnection angle can show us 
whether there is a relationship between the velocity of fan-shaped jets and the reconnecting component 
of magnetic field. For instance, if the reconnection angle is 7r/2, the maximum reconnecting component 
of magnetic field is only Bj„j sin (7r/4) (for other reconnection angle case, Bi„i sin (6'/2)), and the 
velocity of the ordinary reconnection outflow ejected by the magnetic tension force has a similar form 
as Vji sin (n/A) (for other reconnection angle case, va sin {0/2)). We know that the dependence of the 
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X 



Fig. 9 A cartoon for explaining the fan-shaped jets. 

velocity on the resistivity v alue, if it is localized, is not so remarkable in 2D magnetic reconnection sim- 
ulations (Ijiang et al.ll2010l) . However, the driving force for these jets seems different from the magnetic 
tension force. Thus, we have to check such a dependence. The dependence on other parameters will be 
studied in the near future. 

The dependence on the resistivity is very simple. As shown in Figures[T6]and[T7] the variables and 
velocity are constant for different resistivity values. One important thing is that there are two kinds of 
velocity we used in these figures. One is the maximum velocity in the computational box which stands 
for the velocity of the ordinary reconnection jets, the other is the maximum parallel velocity (parallel 
to the initial magnetic field) which indicates the speed of fan-shaped jets. Moveover, the reconnection 
becomes faster when the resistivity value increases. We can see that there is no remarkable dependence 
between simulation results (density, pressure, temperature and velocity) and the initial resistivity values. 
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Fig. 10 Velocity distributions (typical case) along x, y, and z directions at different times. 
The Alfven speed is about 1.58. 
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Time=0.00 Time=8.50 




Fig. 11 The panels show the trajectories (dark lines) of a lagrangian fluid element with the 
magnetic field lines (typical case). The blue dashed box is the diffusion region. 
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Fig. 12 Variation of the cosine values of the angle between the velocity and the magnetic field 
direction (left panel). Right panel shows the forces per unit mass and velocity in z direction. 
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Time=0.00 




Fig. 13 The upper left panel shows the trajectory (dark line) of a lagrangian fluid element 
(typical case), the blue dashed box is the diffusion region and a 2D distribution in the bigger 
dark dashed plane in the upper left panel is shown in the upper right panel. The lower panels 
depict the gas and magnetic pressure distributions of the small dark box shown in the upper 
right panel at the time = 11.50. The fluid element is marked by the white diamond symbol. 
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Fig. 14 Variations of the reconnection rate (upper left panel), electronic field at the x point 
(upper right panel), energy and mass conservations (lower left panel) and energy release rate 
(lower right panel) in the typical case. and Ey are almost zero, so we can not see them in 
the upper right panel. Wt, Wk and Wm in the lower right panel are the total thermal energy, 
total kinetic energy and total magnetic energy, respectively. 
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Fig. 15 Upper panels show the gas pressure and current density distributions in x — y plane, 
which is also the central part of Figured Other panels depict one dimensional plots along the 
white line in the upper left panel at the time = 14.0. A pair of slow mode shocks are shown 
by the current density distribution and the shadows in the ID distributions. The red lines are 
used to measure the angle (a) between the two shocks. 
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Fig. 16 The maximum and the minimum densities, pressures and temperatures as a function 
of the resistivity values. 
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Fig. 17 The left panel shows the maximum velocity indicating the speed of ordinary recon- 
nection jets) and the parallel velocity (parallel to the initial magnetic field which means the 
velocity of fan-shaped jets) as a function of the resistivity values. The right panel depicts the 
reconnection rates presented by velocity. Note that the reconnection angle (9) is tt in this case. 



Figure [18] depicts 3D visualizations for different reconnection angles. The left panels show gas 
pressure distributions with a velocity field and in the z = cross section plane. The right panels also 
show gas pressure distributions but with a magnetic field and a different view point. The velocity along 
the guide field is weaker than that in the typical case. Figures [19] and [20] show the 2D gas pressure 
distributions, which are similar to that in the typical case. Figure [2T]depicts the maximum and minimum 
densities, pressures and temperatures as a function of the reconnection angle. These values are measured 
at the time when the reconnection rate reaches the maximum. If we decrease the reconnection angle, the 
difference between the maximum and the minimum values becomes smaller, and finally the maximum 
value will be equal to the minimum one and being the initial value 1 . Linear relation is good in this 
dependence. The left panel of Figure [22] shows the maximum velocity in the computational box which 
stands for the velocity of the ordinary reconnection jets, the other is the maximum parallel velocity 
(parallel to the initial magnetic field) which indicates the speed of fan-shaped jets. From the left panels 
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in Figures [17] and |22] we found that the speed of fan-shaped jets is around half of the speed of ordinary 
reconnection jets, and it also depends on the reconnection angle (the reconnecting component of the 
magnetic field). The right panel of Figure |22] shows the reconnection rates for different angles. The 
maximum reconnection rate is about 0.06. The smaller angles correspond to a longer evolution time, 
which is reasonable. 

5 DISCUSSION AND SUMMARY 

We perform three dimensional resistive MHD simulations to study the magnetic reconnection using an 
init ially shearin g magnetic field configuration, which is similar to the papers by Pontin et al. (2005i) 
and lUgail(l2010l) . It is shown that there are two types of reconnection jets: the ordinary reconnection jets 
which move along the ambient magnetic field lines, and the fan-shaped jets which move along the guide 
field lines. In this paper we analyse the fan-shaped jets in detail. This may provide a new understanding 
of the dynamic phenomena in astr ophysics. One of the po ssibility is that it can be applied to in terpret 
the solar sunsport penumbral jets (iKatsukawa et al.|[2007b or anemone jets dShibata et al.ll2007l) in the 
solar chromosphere. 

Actually, before the fan-shaped jets, there is a slow mode wave propagating along the magnetic field 
lines. Because of its fast speed this wave can be seen as the front of the fan-shaped jets. But it is very 
hard to distinguish which is the slow mode wave and which are the jets. There is no clear boundary 
between them. The part which has gas pressure increasement but with small velocity is the slow mode 
wave as shown in Figure|2] Of course, the ordinary jets connect the fan-shaped jets, the slow mode wave 
also exists at the front of ordinary reconnection jets. 

The driving force for the fan-shaped jets is another interesting point. As there is no gravity in our 
simulations, the analysis of the Lagrangian fluid elements is easy to understand. Before Lagrangian 
fluid elements move along the magnetic field lines, both the gas pressure gradient force and the Lorentz 
force can drive the elements. The key point is that the magnetic tension force dominates the element 
movement first and then the element is affected by the gas pressure gradient force. In case 1 listed in 
Table [2] with a small total grid points, we incidentally found that a fluid element is only accelerated by 
the Lorentz force. That is due to the initial position of that element being at the upper edge (outflow side) 
of the diffusion region. It goes out of the diffusion region soon, so that only the Lorentz force affects 
this element. 

The dependence of the simulation variabes (namely density, pressure, temperature and velocity) on 
resistivity and reconnection angle is reasonable. The measurement of reconnection rate is simple in our 
paper and needs to be improved. However, from this method, people can get a general understanding 
about the 3D magnetic reconnection rate. It is roughly but useful, at least we can know when the mag- 
netic reconnection reaches the maximum. Another disadvantage in our simulation is the limitation of 
the total grid points. A finer grid mesh or a higher order finite difference scheme or algorithm is under 
consideration. 

If the effect of gravity is considered, the density stratification effect in the solar atmosphere can 
greatly increase velocity amplitude of the upward fan-shaped jets when it propagates up to the low 
density region (like in the solar upper chromosp here and the corona) and eventually becomes a shock 
or jets dShibata & Suematsulll982HShibatalll982l) . The extended development of this study with gravity 
considered is in progress. 

In summary, we give the conclusions as follows: 

1. From 3D MHD simulations we have studied a type of jets which move along the guide magnetic 
field. Because of the rotation of the initial magnetic field these jets have a fan-shaped structure. Of 
course, the ordinary reconnection jets ejected by the magnetic tension force also exist. 

2. The driving forces of these fan-shaped jets are Lorentz force and gas pressure gradient force. The 
magnetic pressure gradient force drives the Lagrangian fluid elements first, and then the force on the 
elements is dominated by the gas pressure gradient in the later stage. 

3. The pressure, density, temperature and velocity do not sensitively depend on the initial resistivity 
value. 
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Fig. 18 3D visualizations of the gas pressure distributions at different times (The reconnection 
angles for these cases are 37r/4, 7r/2 and 7r/3 from top to bottom). The bright colors and high 
opacity stand for the high gas pressure and the low value is transparent. A z = cross section 
and a velocity plane are used in the left panels. Magnetic field and a different view points are 
used in the right panels. 
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Fig. 19 2D distributions of gas pressure in the x — y, y ~ z and x — z planes for different 
cases (The reconnection angles for these cases are 37r/4, tt/2 and 7r/3 from top to bottom, 
respectively). The arrows indicate the magnetic field. 



4. The ratio between the velocity of the fan-shaped jets and the velocity of the ordinary reconnection 
jets is almost a constant being about 0.5 for different reconnection angles. 
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Fig. 20 2D distributions of gas pressure in the y~z plane for different cases (The reconnection 
angle for these cases are 37r/4, tt/2 and 7r/3 from top to bottom, respectively). The arrows 
indicate the velocity field. 
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